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-Abstract- 

We study networks obeying time-dependent min-cost path metrics, and present novel oracles 
for them which provably achieve two unique features: (i) subquadratic preprocessing time and 
space, independent of the metric’s amount of disconcavity; (ii) sublinear query time, in either 
the network size or the actual Dijkstra-Rank of the query at hand. 
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[T] Introduction 

Concurrent technological infrastructures (e.g., road networks, social networks, e-commerce 
platforms, energy-management systems) are typically of very large scale and impose as a 
routine task the computation of min-cost paths in real-time, while their characteristics usu¬ 
ally evolve with time. The large-scale and real-time response challenges have been addressed 
in the last 15 years by means of a new algorithmic trend: the provision of oracles. That 
is, data structures created by appropriately selecting precomputed information (summaries) 
and which subsequently support query algorithms with real-time responses. The quality of 
an oracle is assessed by its preprocessing space and time requirements, the time-complexity 
of the query algorithm and the approximation guarantee (stretch). Numerous oracles have 
been proposed and analyzed (see e.g., [Tl [Ml[Ml[?7l and references therein) for 
large-scale, mostly undirected networks, accompanied by a static arc-cost metric. In tandem 
with oracles, an equally important effort (with similar characteristics) has also emerged in 
the last 15 years under the tag of speedup techniques, for approaches tailored to work ex¬ 
tremely well in real-life instances (see e.g., |3] and references therein). 

The temporality of the network characteristics is often depicted by some kind of pre¬ 
determined dependence of the metric on the actual time that each resource is used (e.g., 
traversal speed in road networks, packet-loss rate in IT networks, arc availability in social 
networks, etc). Perhaps the most typical application scenario (motivating also our work) is 
route planning in road networks, where the travel-time for traversing an arc a = uv (mod¬ 
eling a road segment) depends on the temporal traffic conditions while traversing uv, and 
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thus on the departure time from its tail u. This gives rise to time-varying network mod¬ 
els and to computing min-cost (a.k.a. shortest) paths in such networks. Several variants 
of this problem try to model time-variation of the underlying graph structure and/or the 
arc-cost metric (e.g., dynamic shortest paths, parametric shortest paths, stochastic shortest 
paths, temporal networks, etc). In this work, we consider the case in which the cost vari¬ 
ation of an arc a is determined by a function D[a\, which is considered to be a continuous, 
piecewise linear (pwl|^and periodic function of the time at which the resource is actually 
being used [3 m HI US]- In case of providing route plans in time-dependent road networks, 
arc-costs are considered as arc-travel-time functions, while time-dependent shortest paths as 
minimum-travel-time paths. The goal is then to determine the cost {minimum-travel-time) 
of a shortest path from an origin o to a destination d, as a function of the departure-time to 
from o. Due to the time-dependence of the arc-cost metric, the actual arc-cost value of an 
arc a = uv is unknown until the exact time > to at which uv starts being traversed. 

Problem Setting and Related Work. Two variants of the time-dependent shortest path 
problem have been considered in the literature: TDSP{o,d,to) (resp. TDSP{o,-k,to)) fo¬ 
cuses on the one-to-one (resp. one-to-all) determination of the scalar cost of a minimum- 
travel-time (shortest) path to d (resp. for all d), when departing from the given origin o 
at time to- TDSP{o,d) (resp. TDSP{o,-k)) focuses on the one-to-one (resp., one-to-all) 
succinct representation of the time-dependent minimum-travel-time path function{s) D[o, d] 
from o to c? (resp. towards all reachable d), and all departure-times from o. 

TDSP{o,d,to) has been studied as early as [ 3 . The first work on TDSP{o,d,to) for 
continuous time-axis was m where it was proved that, if waiting-at-nodes is allowed uncon¬ 
ditionally, then TDSP{o,d,to) is solvable in quasilinear time via a time-dependent variant 
of Dijkstra’s algorithm (we call it TDD), which relaxes arcs by computing the arc costs “on 
the fly”, upon settling their tails. A more complete treatment for the continuous case, con¬ 
sidering various limitations in the waiting-times at the nodes of the network, was provided in 
[13]; an algorithm was also given for TDSP{o,d,to), whose complexity cannot be bounded 
by a function of the network topology. An excellent overview of the problem is provided in 
|19j . Among other results, it was proved that for affine arc-cost functions possessing the 
FIFO property (according to which all the arc-cost functions have slopes at least —1), in 
addition to TDD, a time-dependent variant of the label-correcting Bellman-Ford algorithm 
works. Moreover, if waiting-at-nodes is forbidden and the arc-costs do not preserve the FIFO 
property, then subpath-optimality of shortest paths is not necessarily preserved. In such a 
case, many variants of the problem are also NP-hard m- Additionally, when shortest path 
costs are well defined and optimal waiting-times at nodes always exist, a non-FIFO arc with 
unrestricted-waiting-at-tail policy is equivalent to a FIFO arc in which waiting at the tail 
is not beneficial m- For these reasons, we focus here on instances for which the FIFO 
property holds, as indeed is the case with most of past and recent work on TDSP{o,d,to)- 

The complexity of TDSP{o,d) was first questioned in jSj |7] and remained open until 
recently, when it was proved in |12j that, in case of FIFO-abiding pwl arc-cost functions, 
for a single origin-destination pair (o, d) the space complexity for succinctly representing 
D[o, d] is (1 -T K) ■ where n is the number of vertices and K is the total number of 

breakpoints (or legs) of all the arc-cost functions. Note that K can be substituted by the 


Major car navigator vendors provide real-time estimations of travel-time values by periodically sampling 
the average speed of road segments using the cars connected to the service as sampling devices. The 
most customary way to represent this historic traffic data, is to consider the continuous pwl interpolants 
of the sample points as arc-travel-time functions of the corresponding instance. 
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number K* of concavity-spoiling breakpoints (at which the arc-cost slopes increase) of the 
arc-cost functions. Several output-sensitive algorithms for the exact computation of D[o^d] 
have been presented [HIHIIIIIIS], the most efficient being the ones in El HI]. 

Due to the above mentioned hardness of TDSP{o, d), and also since the time-dependent 
arc-costs are typically only (e.g., pwl) approximations of the actual costs, it is quite natural 
to seek for succinct representations of approximations to D[o,d], which aim at trading-off 
accuracy for computational effort. Several one-to-one (1 -I- e)-approximation algorithms for 
TDSP{o,d) appeared recently in the literature EIH^IISIj the most successful being those 
provided in m- The first one-to-all (l-l-e)-approximation algorithm for TDSP{o,-k), called 
bisection (BIS), was given in |l(i] . It is based on bisecting the (common to all functions) axis 
of departure-times from o and considers slightly stricter assumptions than just the FIFO 
property for the arc-cost metric. BIS requires • log^ calls to TDSP{o,*,to)- 

Note that all one-to-one approximation algorithms for TDSP{o,d) El 1121 HE] demand, in 
worst-case, a comparable amount of calls to TDSP{o,-k,to), just for one od-pair. 

Minimum-travel-time oracles for time-dependent networks (TD-oracles henceforth) had 
received no attention until recently m- A TD-oracle is an offline-precomputed data struc¬ 
ture that allows the evaluation of an upper-approximation A[o, d](to) of D[o,d\{to) (due to 
the hardness of computing and storing D[o, dj), for any possible query (o, d, to) G I^x MxK> o 
that may appear in an online fashion. One trivial solution would be to provide a succinct 
representation of A [o, d] for each pair (o, d) G V xV, for the sake of rapid evaluations in the 
future but at the expense of superquadratic space. Another trivial solution would be to ex¬ 
ecute TDD “on-the-fly” for each new query (o, d, to), at the expense of super linear query-time. 
A non-trivial TD-oracle should aim to trade-off smoothly preprocessing requirements with 
query-response times and approximation guarantees. In particular, it should precompute a 
data structure in subquadratic time and space, and also provide a query algorithm which 
evaluates efficiently (i.e., faster than TDD) A[o, d](to), where A[o, d] must have a provably 
good approximation guarantee w.r.t. D[o, dj. Note that there exists important applied work 
{speedup heuristics) for computing time-dependent shortest paths (e.g., El El HHl HZj) ? which 
however provide mainly empirical evidence on the success of the adopted approaches. 

The TD-oracles in m required 0{ffi ^{K*-\-l)) preprocessing space and time, for 
some constant jd G (0,1), and are able to answer queries (under certain conditions) in time 
0{ffi), for some constant S G (0,1). When K* G o{n), the oracles can be fine-tuned so 
as to assure sublinear query-times and subquadratic preprocessing requirements in n. An 
extensive experimental evaluation of those oracles on a real-world road network is provided in 
m , demonstrating their practicality, at the expense, however, of large memory consumption 
due to the linear dependence of the preprocessing space on K* which can be Q,{n). 

The main challenge addressed here is to provide TD-oracles that achieve: (i) subquadratic 
preprocessing requirements, independently of K*\ and (ii) query-times sublinear, not only 
in the worst-case (i.e., in n), but also in the number F[o, (i](to) of settled vertices when 
executing TDD(o, *, to) until d is settled (Dijkstra-Rank). 

Our Contributions. We address positively the aforementioned challenge by providing: 
(i) A novel and remarkably simple algorithm (TRAP) (cf. Section]^ for constructing one-to- 
many (1-P e)-upper-approximations A[o, d] {summaries) of minimum-travel-time functions 
D[o, d], for all “sufficiently distant” destinations d from the origin o. TRAP requires o(n) calls 
to TDSP{o,-k,to), which is independent of the degree of concavity K* . Its novelty is that it 
does not demand the concavity of the unknown function to approximate, (ii) The TRAPONLY 
and FLAT oracles (cf. Section]^ which exploit TRAP and BIS to construct minimum-travel¬ 
time summaries from randomly selected landmarks to all reachable destinations. The pre- 
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preprocessing space/time 

query time 

recursion budget (depth) r 

m 

K* ■ 

^(5+o(l) 

r £ 0(1) 

TRAPONLY 

„2-/3-Hc(1) 

j^d+o(l) 

r « ^ - 1 

a 

FLAT 

„2-/3-Hc(1) 

72*+o(l) 

r « ^ - 1 

a 

HORN 

^2-/3-Ho(l) 

^«r[o,d](4)^+di) 

r « ^ - 1 

a 


M Table 1 Achievements of oracles for TD-instances with period T = n“, for constant a £ (0,1). 
The stretch of all query algorithms is 1 + e • ■ For all oracles, except for the first, we 

assume that /3 ^ 0. 


processed data structures require subquadratic space and time, independently of K*. FLAT 
uses the query algorithms in m- TRAPONLY needs to extend them in order to recover missing 
summaries for local neighborhoods around a landmark. In both cases sublinear query-times 
are achieved, (in) The HORN oracle (cf. Section]^ which organizes a hierarchy of landmarks, 
from many local landmarks possessing summaries only for small neighborhoods of destin¬ 
ations around them, up to a few global landmarks possessing summaries for all reachable 
destinations. HORN’S preprocessing requirements are again subquadratic. We then devise 
and analyze a novel query algorithm (HQA) which exploits this hierarchy, with query-time 
sublinear in the Dijkstra-Rank of the query at hand. 

Except for the choice of landmarks, our algorithms are deterministic. A recent experi¬ 
mental study m demonstrates the excellent performance of our oracles in practice, achieving 
considerable memory savings and query times about three orders of magnitude faster than 
TDD, and more than 70% faster than those in [Tl]. Table summarizes the achievements of 
the TD-oracles presented here and their comparison with the oracles in m- 
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Preliminaries 


Notation and Terminology. For any integer fc > 1, let [fc] = {1, 2,..., A:}. A time- 
dependent network instance (TD-instance henceforth) consists of a directed graph G = {V, A) 
with \V\= n vertices and \A\ = m G 0{n) arcs, where each arc a £ A is accompanied with 
a continuous, pwl arc-cost function D[a] : ]R> o t K>o- We assume that all these functions 
are periodic with period T > 0 and are defined as follows: Vfc £ N,Vt G [0,T), D[a]{kT -F 
t) = (i[a](t), where d[a] : [0,r) —>• (0, M^] is such that limt'lT d[a]{t) = d[a](0), for some 
fixed integer Ma denoting the maximum possible cost ever seen for arc a. Let also M = 
maxagA Ma denote the maximum arc-cost ever seen in the entire network. Since D[a\ is 
periodic, continuous and pwl function, it can be represented succinctly by a sequence of Ka 
breakpoints (i.e., pairs of departure-times and arc-cost values) defining (i[a]. K — 
is the number of breakpoints representing all arc-cost functions, iTmax = maxagA Ka, and 
K* is the number of concavity-spoiling breakpoints (the ones at which the arc-cost function 
slopes increase). Clearly, K* < K, and K* = 0 for concave arc-cost functions. 

To ease the exposition and also for the sake of compliance with terminology in previous 
works (inspired by the primary application scenario of route planning in time-dependent road 
networks), we consider arc-costs as arc-travel-times and time-dependent shortest paths as 
minimum-travel-time paths. This terminology facilitates the following definitions. The arc- 
arrival-time function of a £ A is Arr[a\{t) = t-\- D[a\(t), Vt £ [0,oo). The path-arrival-time 
function of a path p = (oi,..., Cfc) in G (represented as a sequence of arcs) is the composition 
Arr[p]{t) = Arr[afe](Arr[afc_i](- • • (Arr[ai](t)) • • •)) of the arc-arrival-time functions for the 
constituent arcs. The path-travel-time function is then D[p\(t) = Arr[p\(t) — t. 
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For any (o, d) gV xV, Vo,d denotes the set of od-paths. For any p G Vo,x and q G Vx,d, 
s = p • q G Vo,d is the concatenation of p and q at x. The earliest-arrival-time function 
is Arr[o, d]{to) = min^gp^ ^ {Arr[p\{to)}, ytg > 0, while the minimum-travel-time function 
is defined as D[o,d](to) = mmp^'p^^{D[p](to)} = Arr[o,d]{to) — to- For a given query 
{o,d,to), SP[o,d\{to) = {p G Po,d ■ Arr[p]{to) = Arr[o,d\{to)} is the set of earliest-arrival¬ 
time (equivalently, minimum-travel-time) paths. ASP[o^d\{to) is the set of od-paths whose 
travel-time values are (l-|-£)-approximations of the minimum-travel-time among all od-paths. 

When we say that we “grow a TDD ball from {o,to) ” we refer to the execution of TDD 
from o G V at departure time tg G [0,r) for solving TDSP{o,-k,to) (resp. TDSP{o,d,to), 
for a specific destination d). Such a call, which we denote as TDD(o, *, to) (resp. TDD(o, d, to)), 
takes time C>(m-Tnlog(n)[l-T loglog(l-T dfmax)]) = C>(n log(n) log log(iFmax)]), using pre¬ 
decessor search for evaluating continuous pwl functions (cf. [15]). The Dijkstra-RankT[o, d](to) 
is the number of settled vertices up to d, when executing TDD(o,d,to)- 

t/a = uv G A and [ts,tf) C [0,T), we define upper- and lower-bounding travel-time 
metrics: the minimally-congested travel-time D\uv\(ts,tf) := mint^^[t^^tf){D[uv](tu)} and 
the maximally-congested travel-time D[uv]{ts,tf) := maxt^g[t^_tj.){Z3[Mu](t„)}. If [ts,tf) = 
[0, T), we refer to the static free-flow and full-congestion metrics D, D : A ^ [1, M], respect¬ 
ively. Each arc a G A is also equipped with scalars id [a] and D[a] in these static metrics. 
For any arc-cost metric D, diam{G, D) is the diameter (largest possible vertex-to-vertex 
distance) of the graph. For example, diam{G, D) is the free-fiow diameter of G. 

In our TD-instance, we can guarantee that T > diam{G, D). If this is not the case, 
we can take the minimum number c of consecutive copies of each d[a] as a single arc- 
travel-time function d'[a] : [0, cT) >->• IR>o and D'[a]{t-\-kT') = d'[a\(t), V< G [0,T') such that 
T' = cT > diam{G, Df)- In addition, we can also guarantee that T = n°^ for a small constant 
a G (0,1) of our control. If T n“, we scale the travel-time metric by setting D” = ■ D 

(e.g., we change the unit by which we measure time, from milliseconds to seconds, or even 
minutes) and use the period T” = n“, without affecting the structure of the instance at all. 
From now on we consider w.l.o.g. TD-instances with T = n“ > diam{G,D). 

For any v G V, departure-time ty G M>0) integer F G [n] and R > 0, B[v,F]{tv) 
{B[v; i?](t«)) is a ball of size F (of radius R) grown by TDD from {v, ty), in the time-dependent 
metric. Analogously, B\v,F] {B[v, R]) and B[v, F] {B[v,R]) are, respectively, the size-E 
(radius-i?) balls from v in the free-fiow and fully-congested travel-time metrics. 

A pair of continuous, pwl, periodic functions A[o, d] and A[o, d]), with a (hopefully) small 
number of breakpoints, are (1 -T e)-upper-approximation and (1 -T e)-lower-approximation of 
D[o,d\, ifVto > 0, <A[o,d]{to) < D[o,d\{to) < A[o,d\{to) <(! + &) ■ D[o,d]{to) ■ 


Assumptions on the time-dependent arc-cost metric. The directedness and time- 
dependence in the underlying network imply an asymmetric arc-cost metric that also evolves 
with time. To achieve a smooth transition from static and undirected graphs towards time- 
dependent and directed graphs, we need a quantification of the degrees of asymmetry and 
evolution of our metric over time. These are captured via a set of parameters depicting 
the steepness of the minimum-travel-time functions, the ratio of minimum-travel-times in 
opposite directions, and the relation between graph expansion and travel-times. We make 
some assumptions on the values of these parameters which seem quite natural for our main 
application scenario (route planning in road networks), and were verified by an experimental 
analysis (cf. Appendix^. Here we only present a qualitative interpretation of them. It is 
noted that Assumptions |2 . 1 1 and 2.2 were exploited also in the analyses in m- 


► Assumption 2.1 (Bounded Travel-Time Slopes). All the minimum-travel-time slopes are 
bounded in a given interval [—Amin, Amax], for given constants Amin G [0,1) and Amax > 0. 
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► Assumption 2.2 (Bounded Opposite Trips). The ratio of minimum-travel-times in opposite 
directions between two vertices, for any specific departure-time but not necessarily via the 
same path, is upper bounded by a given constant C ^ 1- 

► Assumption 2.3 (Growth of Free-Flow Dijkstra Balls). VF G [n], the free-flow ball B\v\F] 
blows-up by at most a polylogarithmic factor, when expanding its (free-flow) radius up to 
the value of the full-congestion radius within i3[u; F], 

Finally, we need to quantify the correlation between the arc-cost metric and the Dijkstra- 
Rank metric induced by it. For this reason, inspired by the notion of the doubling dimension 
(e.g., [U and references therein), we consider some scalar A > 1 and functions /,g : N i—t 
[l,oo), such that the following hold: y{o,d,to) G V x V x [0,T), (i) F[o,d](fo) < /(n) • 
{D[o, d]{to))^, and (ii) D[o,d\{to) < g{n) ■ (F[o, d]{to)Y^^- This property trivially holds, e.g., 
for A = 1, f{n) = n, and g{n) = max^gA {-D[a]}. Of course, our interest is for the smallest 
possible values of A and at the same time the slowest-growing functions f(n),g(n). Our last 
assumption quantifies the boundedness of this correlation by restricting A, f{n) and g(n). 

► Assumption 2.4. There exist A G f{n),g{n) G polylog(n) s.t. the follow¬ 

ing hold: (i) F[o,d](fo) < /(n) • {D[o,d]{to))^, and (ii) D[o,d]{to) < gin) ■ (F[o,(f](fo))^/^. 
Analogous inequalities hold for the free-flow and the full-congestion metrics I) and D. 

Note that static oracles based on the doubling dimension (e.g., H) demand a constant 
value for the exponent A of the expansion. We relax this by allowing A being expressed as a 
(sufficiently slowly) growing function of n. We also introduce some additional slackness, by 
allowing some divergence from the corresponding powers by polylogarithmic factors. 

In the rest of the paper we consider sparse (m G Oin)) TD-instances, compliant with 
Assumptions |2.1[|2.2[|2.3| and 2.4| For convenience, the notation used throughout the paper 
is summarized in Appendix [A A review of the oracles in m is presented in Appendix [C| 

[J] The TRAP approximation method 

We now introduce the trapezoidal (TRAP) algorithm, a 
novel algorithm for computing one-to-many (1 -F e)- 
upper-approximations A[£, z;] : [0,T) i—>■ ]R>o of min- 
cost functions D[£,v], from a vertex £ towards all 
sufficiently distant destinations (typically, £ will be 
a landmark). TRAP is remarkably simple and works 
as follows. First, [0,r) is split into a number of 
consecutive length-r subintervals, where r is a tun¬ 
ing parameter to be fixed later. Then, for each such 
interval [ts,tf = ts + t) C [0,T), a (1 -F e)-upper- 
approximation of the projection D[£,v] : [ts,tf) i—>■ 

]R>o is computed. Finally, the concatenation of all 
these (1 -F e)-upper-approximations per subinterval 
constitutes the requested (l-Fe)-upper-approximation 
A[£,v] of D[o,d] : [0,T) i-G- M>o. Note that, contrary 
to the BIS approximation algorithm m , no assump¬ 
tion is made on the shapes of the min-cost functions 
to approximate within each subinterval; in particular, 
no assumption is made on them being concave. TRAP only exploits the fact that t is small, 
along with Assumption |2 .1 1 on the boundedness of travel-time slopes. 



H Figure 1 The upper-approximation 
5fc[f, u] (thick orange, upper pwl line), 
and lower-approximation a] (thick 
green, lower pwl line), of the unknown 
function D[l, u] (blue pwl line) within the 
interval F = [ts = {k — 1)t, t/ = fcr). 
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We now describe the upper- and lower-approximations of I?[o, (f| that we construct 
in a subinterval Ik = [ts = (k — l)r, t/ = kr) C [0,T), k G from a vertex 

£ € V towards some destination v € V. The quality of the upper-approximation de¬ 
pends on the value of r and the delay values at the endpoints of Ik, as we shall ex¬ 
plain shortly. TRAP computes the following two functions of D[£,v] (cf. Fig. Vt G 
h, Sk[£,v]{t) = min{ D[£,v]{tf) + A min^/ j A^naK^s Amaxt } and 

—fc [^5 ^](^) — max * 1 ^ 'c](ty) A^axt/ Aj^axt , u](ts) -t- A^ints Aaaint } and con¬ 

siders them as the upper- and lower-approximating functions of D[i,v] within Ik- The 
correctness of this choice is proved in the next lemma, which follows by Assumption |2 .1 [ 

► Lemma 1. Sk[£,v]{t) and Sk[£,v](t) upper- and lower-approximate D[£,v](t) within Ik- 

Let {t^,Dm) {tm,Dm) be the intersections of the legs in the definition of 5^[£,v\ and 
5k[£,v\, respectively. The maximum additive error MAE{Ik) for 5k[£,v\ in Ik (i.e., the 
length of the purple dashed line in Fig. is MAE{Ik) := maxjg/,^ {Sk[£,v]{t) — 5^, [£, u] (t)} = 
~Ek[£, v\{t^)—5^k[^^ ■ The following lemma proves that, for r sufficiently small, MAE{Ik) 

cannot be large. It also provides a sujficient condition for the value of r so that 5k[£,v\ is 
indeed a (1 -P £)-upper-approximation of D[i,v] within Ik- 

► Lemma 2. V(.£,u) G L x V, t/k G [[ 7 ]] and Ik = [(fc — l)r, fcr), the following hold: (1) 

MAE[£,v]{Ik) < Amax • t; (2) 5k[£,v] is a (1 -\- e)-upper-approximation of D[£,v] within Ik, 
if [ D[t,v\{ts) > (Amin + ■ T ] V [ D[£,v\{tf) > (1 -p b) Amax ’ T ] 

For given r > 0 and £ G L, the set of faraway destinations from £ is V[£]{t) = {v G V : 
t[£,v\ > t} . t[£,v] = — is a sufficient r-value for Jk[£,v\ being (1 -P e)-upper- 

approximation of D[£,v\ within Ik = [{k — 1)t[£, u], u]) (cf. Lemma [^. The next 

theorem proves that TRAP provides a (1 -P e)-upper-approximation A[£, u] for all faraway 
destinations from £, and also estimates the preprocessing requirements of the algorithm. 

► Theorem 3. Fix £ G L, F > f{n), and t G (0,T) s-t- |I^[^](t)| = n — F- Let 

T* = min,,g 7 [q(T) I Vu G V[£]{t), A[£,u] is the concatenation of all the 

upper-approximating functions Sk[£,v] that TRAP returns per subinterval Ik = [ ts,. = {k — 
l)r* , tf,, = min{fcr*,r} ) : k G [[:^]] • Then, \/v G V[£]{t), A[£, v] is a {1-\-e)-upper- 
approximation of D[£,v] in [0,r), requiring at most 2 1";^] breakpoints- The number of calls 
to TDSP{£,-k,t) for their construction is < 1 + ^ 0{n°’) . 

Proof of Theorem [31 r* is the appropriate length for the subintervals which assures that 
TRAP returns (l-p£)-upper-approximations for all faraway destinations from £- By definition 
it holds that r* > r. Since F > f{n), it holds that TRAP does not consider destinations at 
free-fiow distance less than 1. To see this, fix any v gV s.t. D\£, v] < 1. By Assumption |2.4[ 
r[^w] A fin) ■ D[£,v]^ < fin) < F - Thus, we can be sure that v ^ tT[^](r). Since 
T = n“, we conclude that ^ ^ We proceed now with the 

analysis of TRAP. [0,T) is split into consecutive length-r* subintervals. Lemma 
assures that for each Ik = [kr*, (A:-P1)t*) an upper-approximating function 5k[£, v] of D[£, u] 
is determined, for each v G V[£]iT)- The concatenation of all these functions constitutes 
the upper-approximating function A[£, u] for D[£,v\ within [0,r). Since t[£,v\ > t* =P 
I^£,v] > (1 + y) AmaxT*, we deduce (cf. Lemma that, for all v G IP[^](r), the produced 
upper-approximations within the consecutive length-r* intervals are (1-P £)-approximations 
of D[£,v]- TRAP preprocesses £ G L (concurrently for all v G V[£]iT)) by making G 

0(n“) calls to TDSPi£, *, t), to sample the endpoints of all the [length-r* subintervals. 
For storing A[£, u], it needs at most 2 breakpoints (there is at most one intermediate 
breakpoint ii-^nDm) per subinterval). ◄ 
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Oracles with fully-informed landmarks 

In this section we describe two novel oracles, with landmarks possessing summaries for all 
reachable destinations, excluding possibly a small neighborhood around them. We start with 
a random landmark set L Cuar(p) V ^ he., we decide independently and uniformly at random 
whether each vertex is a landmark, with probability p = n~‘^ for a constant w S (0,1). We 
consider as faraway vertices from £ £ L, all the vertices at free-flow distance at most R = T^ 
from it, for a constant 6 £ ( 0 , 1 ) to be determined later. F = max^g/, ,R]|} is the 

maximum number of faraway vertices from a landmark. The next lemma shows that the 


main parameters we should consider w.r.t. a TD-instance are A (cf. Assumption 2.41 and 


a £ (0,1) s.t. T = n“. All the other parameters essentially adjust their values to them. 

► Lemma 4. For v £ (0,1) s.t. T = diam{G, , 9 £ (0,1) s.t. | £ 0(1) and X,f,g 


defined as in Assumption 2.4 the following hold: (i) ^ = a±o(l), and (ii) F £ ^ 


The TRAPONLY oracle. A first attempt towards avoiding the dependency of the prepro¬ 
cessing requirements on K* , is to develop an oracle whose preprocessing is based solely on 
TRAP. The preprocessing of this oracle (we call it TRAPONLY) first splits the entire period 
[0,T) into consecutive subintervals of length r = — > 0 each. It then calls TRAP 

for each landmark £ £ L, which guarantees (1-P £)-upper-approximations for all the faraway 
destinations v £ V[£](t) (cf. Theorem]^. As for the faraway destinations from t, we let 
their distances from £ be computed by the query algorithm that we are using (by growing 
a TDD ball from £). In particular, the query algorithm of TRAPONLY is an appropriate vari¬ 
ant of RQA, we call it RQA+, which additionally grows a small TDD ball of size Fpolylog(F) 
(cf. Assumption 2.3) in order to compute the actual travel-times towards their faraway 


destinations. The following theorem analyzes the performance of TRAPONLY. 

► Theorem 5. The expected time o/RQA+ and the preprocessing requirements o/TRAPONLY 


are: E{Qrqa+} G 0 |^„--+max{<..|}+a(i)^ 


Tti 


G C>( 




+o(l)^ . 


Proof of Theorem |5j During the preprocessing, TRAPONLY makes [y] < 1 -P _ 

1 -P T^~^(l -P l/e)Amax calls of TDD(.£,t), for departure-times t £ |0,t,2t, ..., —1} and 

landmarks £ £ L, where the equality comes from Lemma Therefore, the preprocessing¬ 
time is dominated by the aggregate time for all these TDD probes. Taking into account that 
each TDD probe takes time 0(nlog(n) loglog(iLmax)) and that \L\ = pn = landmarks, 
by using Lemma W we get the following: Ptraponly = n}~‘^ ■ • nlog(n) loglog(n) G 


2-bj+: 


41 + 0 ( 1 ) 1 + 


Og log(Tl)+log log log(Tr) 


= 7 j 2 -a;+a (i-e)-i-o(i) ^ Calculations are analogous for 
the required preprocessing space: For all landmarks £ £ L and all their faraway des¬ 
tinations V £ V[£](tY the total number of breakpoints to store is at most S'traponly = 
2 pn^ G n^“‘^++ic4+o(i)]+c(i) _ ^ 2 -a;-i-a'(i-s)+o(i) _ query-time complexity of 

RQA+, recall that the expected number of TDD balls that it grows is (1/p)^ . Additionally, 
RQA+ grows (1/pY TDD balls from the corresponding closest landmarks. Each ball from a 
new center costs 0((1/p)\og(l/p)). Each ball from a landmark costs 0(Fpolylog(F)) G 
Thus, the expected query-time is upper-bounded as follows: E{(5 rqa+} G 
0((l/ pY[( 1/p)^og(l/p) + Fpolylog(F)]loglog(iY,,ax)) = 0 (n™-’^{- 4 i+c(i)]eM) . ^ 

The next corollaries are parameter-tuning examples showcasing the trade-offs among the 
sublinearity of query-time, the subquadratic preprocessing requirements and the stretch. 


► Corollary 6. For 5 £ (a, 1), fi £ (0,a'^iy], uj = — j, 9 = ^ and r = 


i"' 

'S'traponly, Ttraponly G P+dA ^ £ {Qroa+I G fAe stretch is 1 -P £ • ■ 


a -\-0 
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► Corollary 7. For any integer k > 2, let rjik) = 


log(fc/(fc-l)) 

log(l+e/i/i) 


0 , 


r)(fe)+2 


Then TRAPONLY achieves stretch 1 + k ■ e with S- 


— 1, (5 G (0,1) and jS G 

^TRAPONLY G n^ 


TRAPONLY, -f^TRAPONLY 


and EIQroa+I G by scaling the TD-instance so that T = n°‘, for a = ^ • 

The FLAT oracle. Our second attempt, the FLAT oracle, provides preprocessed information 
for all reachable destinations from each landmark. In particular, it uses the query algorithm 
RQA m The preprocessing phase of FLAT is based on a proper combination of BIS and 
TRAP for constructing travel-time summaries. Each landmark i G L possesses summaries 
for all reachable destinations: BIS handles all the (at most F = max^gi i?]|}) nearby 
destinations in whereas TRAP handles all the faraway destinations of V\B[£; R], The 

space requirements for the summaries created by TRAP are exactly the same as in TRAPONLY. 
As for the summaries computed by BIS, we avoid the linear dependence of BIS on K* by 


assuring that F is sufficiently small (but not too small) and exploiting Assumption 2.3 which 
guarantees that the involved subgraph B_'[i; E] in the preprocessing phase of BIS on behalf 
of i has size 0{F polylog(F)). The next lemma shows exactly that BIS is only affected by 
the concavity-spoiling breakpoints of arc-travel-time functions in B^[i;F], rather than the 
entire graph. 

► Lemma 8. V(£,t>) G Lx B[i;F]yu G V\B'[i;F]yt G [0,T), D[i,v]{t) < D[i,u](t). 

Proof of Lemma |8l From the definitions of the involved free-flow and full-congestion Dijk- 
stra balls, the following holds: D[£,v]{t) < D[£,v] < R[£] < D[£,u] < D[£,u]{t). ◄ 

The following theorem summarizes the complexities of the FLAT oracle. 

► Theorem 9. The query-time Qrqa and the preprocessing time Pflat o,nd space S'flat o/FLAT 

are: EIQr^a} G C>(n‘^(’'+i)+o(i)) and Pflat , S'flat G . 

Proof of Theorem [51 BIS requires space at most polylog(E), since by Lemma the 
involved graph only contains F polylog(A) vertices and concavity-spoiling breakpoints at 
the arc-travel-time functions. For the faraway vertices of F \ Bl£-,F], by setting r = 
(i+i/'f)A —) trap provides (1 -F £:)-approximate summaries, since Vu G F \ BIT, i?] the suffi¬ 
cient condition of Theorem]^ holds: D[£, v] > i? = (1 + 1/e) Amax^ • Thus, we conclude that 


Sflat G pn 


F^ polylog(F) -F ^(i+iA)Amaxn 


-F„i+«-(i-e)[i+o(i)]] 


= f^i-“+[i+o(i)] max{ 20 / 1 ^ , i-i-a(i-e) }+o(i) ^ f{n),g{n) G polylog(n). 

The next corollaries are parameter-tuning examples to showcase the effectiveness of FLAT. 
► Corollary 10. If 5 G (a,l), /3 G (o, a; = 4 ,, r = J " 


and 9 = then FLAT has Tflat,Sflat G h+dB ^ E{(5rqa} G and stretch 

(l+e/PY+^ 


1 + £ ■ 




Corollary 11. For integer k > 2, letri{k) = 


log(fc/(fc-l)) 

logil+e/V") 


— 1 and 5 G (0,1). FLAT achieves 


a target stretch 1 -F fc • e with preprocessing requirements and expected query-time 

_ 25 

F(fe)+2]-(2- 


71 > 5 + 0 ( 1 )^ hy scaling the TD-instance so that T = n°‘ for a = \„(y\+2U2-v)-Si^ ’ i 0- 


Comparison 0 /TRAPONLY and FLAT. Both TRAPONLY and FLAT depend on the travel-time 
metric, but are independent of the degree of disconcavity K* . On one hand, TRAPONLY is a 
simpler oracle, at least w.r.t. its preprocessing phase. On the other hand, FLAT achieves a 
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better approximation for the same TD-instance and anticipations for sublinear query-time 
and subquadratic preprocessing requirements . This is because, as /3 0, FLAT 

guarantees a recursion budget r of (roughly) ^ — 1, whereas TRAPONLY achieves about half 
this value and r has an exponential effect on the stretch that the query algorithms achieve. 


The HORN oracle 


We now describe and analyze the Hierarchical OR- 
acle for time-dependent Networks (HORN), whose 
query algorithm is highly competitive against 
TDD, not only for long-range queries (i.e., having 
Dijkstra-Rank proportional to the network size) 
but also for medium- and short-range queries, 
while ensuring subquadratic preprocessing space 
and time. The main idea of HORN is to preprocess: 
many landmarks, each possessing summaries for a 
few destinations around them, so that all short- 
range queries can be answered using only these 
landmarks; fewer landmarks possessing summar¬ 
ies for more (but still not all) destinations around 
them, so that medium-range queries be answered 
by them; and so on, up to only a few landmarks 
(those required by FLAT) possessing summaries for 
all reachable destinations. The area of coverage 
C[P\ C V of £ is the set of its nearby vertices, 
for which £ possesses summaries. £ is called in¬ 
formed for each v € C[£], and uninformed for each 
V G V \C\£]. The landmarks are organized in a 
hierarchy, according to the sizes of their areas of 
coverage. Each level Li in the hierarchy is accom- 


\ 



Rt?s(G3 




H Figure 2 Demonstration of execution 
of HQA. Dashed circles indicate areas of 
coverage. Solid circular stripes indicate 
the rings of the corresponding levels in the 
hierarchy. Landmark li^o is uninformed 
and although informed, comes too 

early. ^ 2,0 is both informed and within the 
ring of its own level, leading HQA to deduce 
that the appropriate level is i = 2. 


panied with a targeted Dijkstra-Rank Ni G [n] , and 

the goal of HORN is to assure that Li should suffice for RQA to successfully address queries 
{o,d,to) with r[o, (i](to) < Ni, in time o{Ni). The difficulty of this approach lies in the 
analysis of the query algorithm. We want to execute a variant of RQA which, based on a 
minimal subset of landmarks, would guarantee a (1 -F tT(r))-approximate solution for any 
query {o,d,to) (as in FLAT), but also time-complexity sublinear in r[o, d](to)- We propose 
the Hierarchical Query Algorithm (HQA) which grows an initial ball from (o^to) that stops 
only when it settles an informed landmark £ w.r.t. d which is at the “right distance” from o, 
given the density of landmarks belonging to the same level with £. HQA essentially “guesses” 
as appropriate level-f in the hierarchy the level that contains £, and continues with the 
execution of RQA with landmarks having coverage at least equal to that of £ (cf. Fig. [^. 


Description of HORN. We use the following parameters for the hierarchical construction: 
(i) k G 0(loglog(n)) determines the number of levels (minus one) comprising the hierarchy 
of landmarks, (ii) 7 > 1 determines the actual values of the targeted Dijkstra-Ranks, one 
per level of the hierarchy. In particular, as 7 gets closer to 1, the targeted Dijkstra-Ranks 
accumulate closer to small- and medium-rank queries, (iii) (5 G (0,1) is the parameter 
that quantifies the sublinearity of the query algorithm (HQA), in each level of the hierarchy, 
compared to the targeted Dijkstra-Rank of this level. In particular, if Ni is the targeted 
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Dijkstra-Rank corresponding to level-i in the hierarchy, then HQA should be executed in time 
0[{NiY), if only the landmarks in this level (or in higher levels) are allowed to be used. 

Preprocessing of HORN. Vi S [k], we set the targeted Dijkstra-Rank for level-i to Ni = 
jiii -i)/7 ^ Then, we construct a randomly chosen level-i landmark set Li Cuar(pi) ^where 
Pi = Ni -i)/[(’'+i)7 Each ii G Li acquires summaries for all (and only 

those) V G C\ii], where C\ii\ is the smallest free-flow ball centered at £i containing Ci = 
Ni-n^' = +4* vertices, for a sufficiently small constant ^i > 0. The summaries to the 

Fi = cf nearby vertices around £i are constructed with BIS; the summaries to the remaining 


_ e _ l+a 

2-\-Oiy 


1 2 
2 ’ 2 + 1 / 


IS 


Ci — Fi faraway vertices of £i are constructed with TRAP, where X = ^ = 
an appropriate value determined to assure the correctness of FLAT w.r.t. the level-i of the 
hierarchy. An ultimate level Lk+i Cuar(pfc+i) ^ of landmarks, with pk+i = assures 

that HORN is also competitive against queries with Dijkstra-Rank greater than . 

We choose in this case Ck+i = Nk+i = n, Fk+i = and C[£k+i] = V, V£k+i G Lk+i- 

Description o/HQA. A TDD ball from (o, to) is grown until d is settled, or the (AS'(7)-criterion 
or the (AL+-criterion is fulfilled (whichever occurs first): 

o Early Stopping Criterion (ESC): £o G L = is settled, which is informed 

{d G C[£o]) and, for (/?>!, ^[A,dKt;,+D[oA](tD) > g.) . i) ^ ^ 

o Appropriate Level of Hierarchy (ALH): For some level i G [k] of the hierarchy, the 
first landmark £i^o & Li is settled such that: (i) d G C[£i^o\ {£i,o is “informed”); and (ii) 

< T[o,£i^o\{to) < ln(n) • {£i^o is at the “right distance”). In that case, HQA 


concludes that i is the “appropriate level” of the hierarchy to consider. Observe that the 
level-(fc + 1) landmarks are always informed. Thus, if no level-(< k) informed landmark is 
discovered at the right distance, then the first level-(fc -F 1) landmark that will be found at 
distance larger than ln(n) • will be considered to be at the right distance, and then 

HQA concludes that the appropriate level is k + 1. 

If d is settled, an exact solution is returned. If (ESC) causes termination of HQA, the value 
D[o,£o]{'to) + ^[£o, d]{to + D[o, £o]{to)) is reported. Otherwise, HQA stops the initial ball due 
to the (ALH)-criterion, considering i G [fc -F 1] as the appropriate level, and then continues 
executing the variant of RQA, call it RQAj, which uses as its landmark set Mi = 

Observe that RQAj may fail constructing approximate solutions via certain landmarks in Mi 
that it settles, since they may not be informed about d. Eventually, HQA returns the best 
od-path (w.r.t. the approximate travel-times) among the ones discovered by RQAj via all 
settled and informed landmarks £. Theorem [T^ summarizes the performance of HORN. 


► Theorem 12. 


Consider any TD-instance with X G o 

e-(r-l-l) 


(V logToSn) ) gin) Jin) G 

polylog(n) (cf Assumption 2.4'l- For p = and k G C>(loglog(n)), let 

C* G ( (1 + A) • loglog(n) -F Alog (^1 -F iny—) /log(n) , 1 - 7 "*), for all i G [fc]. Then, 
for any query (o,d,to) s.t. Ai»_i < r[o, (i](to) < Aj* for some i* G [fc -F 1], any 5 G 
(a, 1), 13 > 0, and r = ^ achieves E{Qhqa} G {Ni.)^+<'^\ 

Thorn , •S'horn G and stretch 1 -F £ , with probability at least 1 — 0(+. 
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I A I Summary of Notation 


Symbol 

Description 

[n] 

The set of integers {1, 2,..., n}, for any n G N — {0}. 

G = (y,A) 

The graph representing the underlying structure of the road network, 
n = \V\ and m = |kl|. 

diam{G, D) 

The diameter of G under an arc-cost metric D. 

'Po,d 

Set of od-paths in G. 

p»q 

The concatenation of the ua-path p with the sn-path q at vertex x. 

ASP[o,d\(to) 

Set of (1 + e)-approximations of minimum-travel-time od-paths in G, for given 
departure-time to > 0. 

5P[o, d\(to) 

Set of minimum-travel-time od-paths in G, for given departure-time to > 0. 

B[v\{U) 

A ball growing from (v,tv) E V x [0,T), in the time-dependent metric, until 
either the destination d or the closest landmark £v from {v,tv) is settled. 

B[v,F\{U) 

A ball growing from {v,tv) E V x [0,T), in the time-dependent metric, of size 

F e N. 

B[v, F] / Blv, F] 

A ball growing from v E V, m the full-congestion / free flow metric, of (integer) 
size F G N. 

B[v,R] / B\v,R] 

A ball growing from v EV ,\a the full-congestion / free flow metric, of (scalar) 
radius R > 0. 

B'[v, F\{U) 

A ball growing from {v,tv) E V x [0,T), in the time-dependent metric, of size 
Fpolylog(F), according to Assumption 2.3 

d-Hii) 

The limited-window arc-travel-time function tor arc a E A, with departure-time 
t E [0,r) for some constant time-period F > 0 (e.g., a single day). 

Ma 

Maximum possible travel-time ever seen at arc a. 

M 

Maximum arc-travel-time ever seen in any arc. 

D\a\{t) 

Periodic arc-travel-time function for arc a E A, with domain t E [0, oo). 

Arr[a\(t) 

The arc-arrival-time function for arc a E A. 

D\o, d\ 

Minimum-travel-time function, from o to d. 

rib;d] 

Dijkstra-Ranks function, from o to d. 

-^max [o, d] / Drain[o,(i\ 

The maximum and minimum value of D[o, d]. 

D[a\ / D[a\ 

Travel-times of a in full-congestion and free-flow metrics, respectively. 

^[o, d] / A[o,d] 

An upper-approximating / lower-approximating function to D[o,d]. 

Arr\o, d\ 

Earliest-arrival-time function, from o to d. 

iu {^v) 

Departure-time from the tail u (arrival-time at the head v) for uv E A. 

TDSP{o, d, to) 

The problem of finding a min-cost od-path, given a departure-time to- 

TDSP{o,*,to) 

The problem of finding a min-cost paths tree, for given departure-time to- 

TDSP{o, d) 

The problem of constructing a succinct representation of min-cost od-paths 
function. 

Ka 

Number of breakpoints in the arc-travel-time function D[a]. 

K 

Total number of breakpoints in the arc-travel-time functions. 

-^max 

The maximum number of breakpoints, among the arc-travel-time functions. 

F* 

Total number of concavity-spoiling breakpoints (i.e., points at which the slope 
increases) in the arc-travel-time functions. 

Ama.x 

Maximum slope among minimum-travel-time functions. 

■^min 

Absolute value of minimum slope among minimum-travel-time functions. 

C 

ratio of minimum-travel-times in opposite directions between two vertices for 
any specific departure-time. 

r 

The recursion budget for RQA and HQA. 

BIS 

The bisection approximation method for minimum-travel-time functions. 

TRAP 

The trapezoidal approximation method for minimum-travel-time functions. 

FCA 

The Forward Constant Approximation query algorithm. 

RQA 

The Recursive Query Algorithm, based on landmarks possessing information 
towards all possible destinations. 

HQA 

The Hierarchical Query Algorithm, based on a hierarchy of landmarks. 

FLAT 

The oracle that uses landmarks possessing summaries towards all possible des¬ 
tinations, and the RQA query algorithm. 

HORN 

The oracle that uses a hierarchy of landmarks with their own subset of destin¬ 
ation vertices, and the HQA query algorithm. 
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B 


Assumptions on the travel-time metric 


In this section, we make a few assumptions on the kind of minimum-travel-time functions in 
the network. All assumptions are quite natural and justified in several application scenarios, 
such as the urban-trafRc road networks, which have motivated this work. Technically, they 
allow a smooth transition from static metrics on undirected graphs towards time-dependent 
metrics on directed graphs. 

The first assumption, called Bounded Travel-Time Slopes, asserts that the partial deriv¬ 
atives of the minimum-travel-time functions between any pair of origin-destination vertices 
are bounded in a known fixed interval: 


2.1 


Assumption 

y X F, VO < ti < t2, 


There exist eonstants Amin G [Oil) ^max > 0 s.t.: V(o,d) G 

Dlo,d](ti)-Dlo,d](t2) ^ A 1 

^ 1 ■'''mini ■'AmaxJ ■ 


The lower-bound of —1 in the minimum-travel-time function slopes is indeed a direct con¬ 
sequence of the FIFO property, which is typically assumed to hold in several time-dependent 
networks, such as road networks. Amax represents the maximum possible rate of change of 
minimum-travel-times in the network, which only makes sense to be bounded (in partic¬ 
ular, independent of the network size) in realistic instances such as the ones representing 
urban-traffic time-dependent road networks. 

The second assumption, called Bounded Opposite Trips, asserts that for any given depar¬ 
ture time, the minimum-travel-time from o to d is not more than a constant C ^ 1 times the 
minimum-travel-time in the opposite direction (but not necessarily along the reverse path). 


Assumption 


2.2 


3C > 1, V(o,d) G y X y, Vt G [0,T), D[o,d\{t) < ( ■ D[d,o]{t). 


This is quite natural in road networks, (i.e., it is most unlikely that a trip in one direction 
is more than 10 times longer than the trip in the opposite direction during the same time 
period). 

A third assumption concerns the relation of the Dijkstra ranks of cocentric balls in the 
network, with respect to the (static) free-flow metric implied by the time-dependent instance 
at hand. Its purpose is to bridge the gap between expansion of graph distances (densities of 
Dijkstra balls) and travel-times in the network. It essentially asserts that, given a particular 
origin-vertex, if one considers a free-flow ball of a certain size (and travel-time radius) and 
then expands further this ball to a larger radius (equal to the full-congestion radius in the 
free-flow ball) then the ball size that we get changes by at most a poly logarithmic factor. 


Assumption [273[ For any vertex £ G V, and a positive integer F, consider the (static) Dijk¬ 
stra ball B[£-, F] around £ under the free-flow metric. Let B\£] = -aiB.y.{D\£, v] : v G B[£-, F]} 
and R[£] = max{D[£,v] : v G B\£', F]} be the largest free-flow and full-congestion travel-times 
from £ to any other vertex in i?[.^;F]. Finally, let lfl[£;F] = {v gV ■. D\£,v]{(),T) < R[£]\ 
be the free-flow ball around £ with the (larger) radius R[£]. Then it holds that |F^[£;F]| G 
0{F ■ polylog(F)). 


The aforementioned assumptions were verified through an experimental analysis on two 
real-world road networks, one concerning the urban-area of the city of Berlin and the other 
concerning the national road network of Germany. Our experimental analysis, presented 
in IISli shows that for the Berlin data set, Amax < 0.062, ( < 1.2, and the maximum ball 
size expansion factor < 6.7, while for the Germany data set Amax < 0.22, ^ < 1.1, and the 
maximum ball size expansion factor <8.3. 

Finally, we need a systematic way to correlate the arc-cost metric (travel-times) with 
the Dijkstra-Rank metric induced by it. For this reason, inspired by the notion of the 
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doubling dimension (e.g., and references therein), we consider some scalar A > 1 and 
functions f,g : N <—>■ [l,oo), such that the following hold: y{o,d,to) G V x V x [0,T), (i) 
T[o,d]{to) < f{n) ■ {D[o,d]{to))^, and (ii) D[o,d]{to) < g{n) ■ (r[o, This property 

trivially holds, e.g., for A = 1, /(n) = n, and g{n) = diam{G,D). Of course, our interest 
is for the smallest possible values of A and at the same time the slowest-growing functions 
/(n),g(n). Our last assumption exactly quantifies the boundedness of this correlation by 
restricting A, f{n) and g{n). 

Assumption [27^ For the graph G = {V,A) and the time-dependent arc-cost metric D that 
we consider, it holds that there exist A G f{n),g{n) G polylog(n) such that 

(i) T[o,d]{to) < f{n) ■ {D[o,d](to))^, and (ii) D[o,d](to) < g{n) ■ (r[o, d](to))^/'^. 

Note that static oracles related to the notion of doubling dimension (e.g., 0 ), demand a 
constant value for the exponent A of the expansion, whereas we allow even a (not too fast) 
growing function of the network size n. The notion of expansion that we consider introduces 
some additional slackness, by allowing some divergence from the corresponding powers by 
poly logarithmic factors of the network size. 


I C I Review of results in [16J 


Review of the approach in flhf - The TD-oracle in [T5] starts by first determining, 
for some p G (0,1), a set L of pn independently and uniformly at random selected land¬ 
marks (vertices acting as reference points). During the preprocessing phase, all (1 -I- e)- 
upper-approximating functions {travel-time summaries) A[£, u] are constructed from each 
landmark i G L towards every reachable destination v G V, using the BIS approximation 
algorithm that keeps bisecting the common axis of departure-times from £, until the desired 
approximation guarantee is achieved in each subinterval, for all destinations. It is proved in 
m that BIS requires 


O 


K* r 

-max < log 

e dev [ 


f T ■ (A„,ax + 1) \ 

\emmt^{D[o,d]{to)} J 


max < log 
dev I 


/ ma.xt^{D[o,d]{to)} \ 1 A 

\mmt^{D[o,d]{to)} ) j J 


calls to TDSP{o,-k,to), for a given origin o gV and all reachable destinations from it. The 
following lemma clarifies this under the lens of our TD-instnace: 

► Lemma 13. BIS requires calls to TDSP{o,*,to) to provide all the sum¬ 

maries of minimum-travel-time functions from a given origin o gV towards all destinations 
at distance at least 1 from o. 

Proof of Lemma 1131 The crucial observation, which is a direct consequence of Assump¬ 
tion [2]^ is that: 

D[o,d\-mo,d\<A„n..-T^ +1 


Exploiting the facts that T = n“, D[o,d] > 1 and A^ax € C>(1), we conclude that BIS 
requires ■ log^ calls to TDSP{o,*,to). ◄ 

Two query algorithms were proposed, FCA and RQA, which provide constant and (1 -I- a)- 
approximations (for constant a > e) to minimum-travel-times, respectively. 

FCA is a simple sublinear-time algorithm for evaluating A[o, d\{to), guaranteeing a con¬ 
stant approximation w.r.t. D[o,d]{to)- In particular, it grows a TDD ball B[o\{to) = 
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{x €V ■. D[o,x\{to) < D[o,io]{to)} from (o,to), until either d or the closest landmark ^o € 
argmmi^L{D[o, i]{to)} is settled. FCA then returns either the exact travel-time value, or the 
approximate travel-time value via £o, A[o,d](to) = D[o,£o]{to) + A[to,d](to + D[o,£o]{to)), 
which is a guaranteed (l-|-£-|-'!/;)-approximation; ijj is a constant depending on e, ^ and Amax, 
but not on the size of the network. RQA improves the approximation guarantee provided by 
FCA, by exploiting carefully a number of recursive accesses to the preprocessed informa¬ 
tion, each of which produces (via calls to FCA) additional candidate od-paths. The tuning 
parameter r G N - the recursion budget - is the depth of the produced recursion tree. 
RQA works as follows: As long as the destination vertex has not yet been discovered in 
the explored area around the origin, and there is still some remaining recursion budget, 
it “guesses” (by exhaustively searching for it) the next vertex Wk at the boundary of the 
current ball, along the (unknown) shortest od-path. Then, it grows a new TDD ball from the 

new center { Wk , 4 = to + D[o,wi]{to) + D[wi,W 2 ]{ti) H-h D[wk-i,Wk]{tk-i) ) until it 

reaches the closest landmark £k to it, at distance Rk = D[wk,£k\{tk)- £k offers an altern¬ 
ative od-path SOLk = Po,wi • • • • • Pwk-i,wk • Qk • Hfc by a new application of FCA, where 
Pwi,wi +1 G SP[wi,Wi+i]{ti), Qk G SP[wk,£k]{tk), and 11^ e ASP[£k,d]{tk + Rk) is the ap¬ 
proximate suffix subpath provided by the oracle. Observe that SOLk uses a longer (optimal, 
if all centers lie on the unknown min-cost path) prefix-subpath Po,u,i • • • • • Pwk-i,wk which 
is then completed with a shorter approximate sufhx-subpath Qk • . It is proved in m 

that the minimum-travel-time over all the discovered approximate od-paths discovered by 
RQA, is a (1 -F cr)—approximation of D[o, d](to), for any constant a > e. The next theorem is 
a consequence of the analysis in m- 

► Theorem 14 ([16]). If a TD-instance with m G 0(n) and compliant with Assumptions 
\2.1\ and \2.2\ is preprocessed using BIS for constructing travel-time summaries from pn land¬ 
marks chosen uniformly-at-random, then the expected values of preprocessing space ^bis 
and time Pbis, cind query time Qrqa for RQA, are: IE{5'bis} G 0{pn^{K* -\- 1)), E{Pbis} G 
0{pn^{K* + l)\og{n)\og\og{K^^Ql), and E {Qrqa} G o(^{l/pf^^ \og{l/p)\og\og{K^^yf}^, 
where r G N the recursion depth in RQA (for r = 0 we get FCA). For the approximation 
guarantees the following hold: FCA returns either an exact od-path, or an approximate od- 
path via a landmark £o s.t. D[o, d](to) < Ro-\-A[£o, d]{to-\-Ro) S (1 + s) ■.D[o, d]{to)-\-tj}-Ro < 
{1 -\- e -\- Ip) ■ D[o,d]{to ), where Ro = D[o,£o]{to) is the minimum-travel-time to the closest 
landmark, and '0 = 1 - 1 - Aniax(l + £)(1 -F 2 ^ -F A^axC) + (1 + £)C ® cost-metric dependent 

constant. RQA returns, for given recursion budget r G N, an od-path that guarantees stretch 
1 -F cr, where a = a{r) < (t|e/^r+i-i ' 

When K* G o(n) these TD-oracles of [121 achieve both sublinear query times and sub¬ 
quadratic preprocessing requirements. Unfortunately, experimental evidence m has demon¬ 
strated that it may be the case that K* G U(n). 


D 


The TRAP approximation method 


In this section, we provide the missing proofs of Section 


Proof of Lemma [T] By Assumption|2.1[ for any departure-time t £ Ik = [ts = {k— 1)t, ty 
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kr) from i and any destination vertex v G V, the following inequalities hold: 


A ^ D[e,v]{t) - D[£,v]{t,) ^ ^ 

— Amin S -^-:- S A 


t - 


_ -^^max 


-Ai„in(t-ts) +D[£,v]{ts) < D[£,v]{t) < D[£,v]{ts) +Amax(t-ts) 


A ^ D[£,v]{tf) - D[£,v]{t) ^ , 

_ , , _ -^^max 

tf -1 

Amin(^/ -t) + D[£,v]{tf) > D[£,v]{t) > -Amax(i/ -t)+ D[£,v]{tf) 


Combining the two inequalities we get the following bounds: Vu € V,Vt € Ik- 

j “Amint + Amints + 'y](ts) 1 < nr/i „,1 ('j-'i < / A^axt — Amax^s + 7? [-^j'f] (tg) 

\A^,^t-A^,^tf + D[£,v]{tf)i- \ -A^,^t + A^intf + D[£,v]{tf) 

Exploiting the fact that each minimum-travel-time function from £ to any destination v G V 
and departure time from Ik respects the above mentioned upper and lower bounds, one 
could use a simple continuous, pwl approximation of D[£,v\ within this interval: 


Vt G 4, 5k[£,v\{t) = min 


-^max^ -I-£)[£,u](ts) - Amaxtg, 1 

+ D[£,v]{tf) + A 

min^/ J 


I.e., we consider the lower-envelope of the lines passing via the point {tg, D[£,v]{tg)) with 
the maximum slope Amax) and the point (tf,D[£,v](tf)) with the minimum slope —Anii„. 
Analogously, we construct a lower-bounding approximation of D[£,v] within Ik' 


Wt G Ik, 5f.[£,v]{t) = max 


^max^ + D[£,v]{tf) - Amaxt/) 1 
+ D[£,v]{ts) + A 

niin^s J 


Figure [^shows the (upper and lower) approximations with respect to D[£,v\ within 


Proof of Lemma |5J Since is the intersection of two lines, it is easy to show that: 

D[£,v]{ts) - D[£,v]{tf) A mints “t- Aaaaxt/ 


t™ = 


= 


A. 


mm ~ -‘^max 


A, 


A. 


mm ~ -‘^max 


A, 


AmaxAl[£, -f AminAl)'^,'c](tj) A 

min ’ -^n 


■^min H” -^max 


■^min H” -^max 


(^/ ^s) 


Analogously, is also the intersection of two lines. Therefore: 

D[£,v]{tf) - D[i,v]{ts) A mint/ T Ajnaxtfi 


tm — 


A 


mm — -ti-max 


A. 


A 


mm — -ti-max 


A. 


j-j AiqiaxT^['t, r’](t/) + Aiqiin-t^[-t, t?] (tg) AnqinA^qax /, , \ 

T'm, = -:-^:-h :- (tf — ts) 


Amin T Amax 


Amin T Amax 
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We start with the upper bound on the maximum absolute error: 


MAE%v]{h)<Dm-D^ 




^min-^^max 


-^min H” -^max -^min H” -^n 

-^max'^['^5 '^] {^s) H" (^/) min-^n 


^min-^^max 


{tf - ts) 
(tf-ts) 


t^min “t” Ajain “1“ A^_ 

(Amax - A^iri)[D[£,v]{tf) - D[£,v]{ts)] + 2K^i^Kmax{tf - tg) 


Aa 


■A, 


/. AbITTIv 
< 


(Amax - Aniin)[A>[^,^](t/) - D[(., v]{t g)] / {t f - tg) + 2AniinA„ 
Amin + A 

(Amax Amin) Amax “1“ 2AminA 


-{tf tg) 


mm ~ -^^max 

max 


Amin “t“ Amax 
Recall now about u] that: Vt € Ik, 


{tf tg^ — Amax'^ 


Sk[i,v]{t) < Sk[£,v]{t) + MAE[e,v]{Ik) < S^.[£,v]{t) + A„ 


< D[£,v]{t)- [1 + 


A„ 


Our goal is to assure that this last upper bound of Sk[£,v]{t) is in turn upper-bounded by 
(1 -I- e) • D[£,v]{t). Based on the expression of 5^[£,v\{t), and exploiting also the fact that 
T > max{t — tg,tf — t}, a sujjicient condition for this to hold, is the following: 

D[i,v]{tg)>{A^in+^)T V D[£,v]{tf) > ^)t 

This sufficient condition is independent of the actual departure time t G Ik, and only depends 
on the travel-time values at the endpoints tg and tf, and also on the length r of the departure- 
times subinterval that we choose. ◄ 

Pe] The FLAT and TRAPONLY oracles 

We start by providing the missing proof of Lemma 

Proof of Lemma |4j All the properties mentioned in the statement of the Lemma, are con¬ 
sequences of Assumption |2.4| We also exploit the fact that the free-flow diameter diam{G, D) 
of G corresponds to the maximum possible Dijkstra-Rank, which is equal to n, assuming 
that the graph is strongly connected, if we use as root of the Dijkstra tree the origin of a 
longest minimum-travel-time path in {G,D). In particular, we proceed with the explanation 
of each property separately: 

(i) Recall that f{n),g{n) G polylog(n) = n^( ) and A G o^^ ■ We start by 

providing the upper bound of T: 

diam{G,D) < g(n)-n^/^ 

=> T = {diam{G,D)Y/’^ <{g{n)f/'' 


= + i[i+A.o(L^ 


£ log log(Tl) 


)] 1 +0(1) 
J] = n i'A 


As for the lower bound of T, we have: 
diam{G,D) > 


T = {dzam{G,D)y^‘'>{j^) 
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(ii) For the upper bound on F we have: 

F = max£gL{|B[^,i?]|} 

where the last step is because we consider instances with A € i s-nd moreover 

due to the fact that we set our (yet unspecified) tuning parameter 0 so that | G 0(1)- 
For the lower bound on F we have: 

E < ff(n) ■ F^/^ 

- \9(n)) (g{.n)y — 

where the last equality is, again, valid since A S | S 0(1). 


E.l Analysis of the TRAPONLY oracle. 

Recall that the preprocessing of the TRAPONLY oracle is based solely on TRAP for computing 
travel-time summaries, while the query algorithm is an appropriate variant of RQA (we call 
it RQA+) which additionally grows a small TDD ball as soon as it settles a new landmark, 
in order to compute “on the fly” the exact minimum-travel-times (rather than evaluating 
preprocessed summaries, which do not exist) towards the nearby destinations from it. The¬ 
orem provides the performance of the TRAPONLY oracle. 

Corollaries and explore the conditions under which sublinear query-time and/or 
subquadratic preprocessing complexities can be guaranteed. We provide here their proofs. 


Proof of Corollary We consider scaled TD-instances with ^ = a for some constant 
a G (0,1). We start with the sublinearity of the query time. For arbitrary constant S G (0,1) 
such that S = max {w • (r -P 1), wr -P ® } , if we choose w = ® then we have <5 = w • (r -P 1) 
w = ^ = 0 = and the sublinearity of the expected query time is guaranteed. 

We continue with the demand for subquadratic preprocessing requirements. Fix now 
some /3 G (0,1). In order to assure preprocessing time and space (roughly) it is 

necessary to assure that 


/3<a; — a-(l — 0) 


S 

r -P 1 


— a -P 


aSiy 

r -P 1 


<tP> 


r < 


5-{l+av) 

a+j3 




► Remark. Note that the query-time performance of RQA+ is equal to that of RQA (cf. [TR]'). 
if w = f. Moreover, P < a^v implies that f. 


Proof of Corollary 1^ We know that, for those instances for which it works, TRAPONLY 
achieves stretch 1 -P e • 


(l+e/^Y+'^-l ■ 


We must assure then that 


(l + ^/^)”+i-l 


< fc =P r > 


log(T^) 

log(l-Pe/V') 


1 = 77(fc) 
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(5- (1 + Qt 


-J — 1 (cf Corollary 


Recall that the maximum recursion budget of TRAPONLY is 
sufficient condition for guaranteeing the required budget for having stretch 1 
the following: 


■ e IS 


6 |. A 
hus 


<5 • (1 -P av) 

a + P 


— 2 > r](k) 


a < 


S-Mk)+2]-l3 


— 7j{k)-\-2 — 6u 




E.2 Analysis of the FLAT oracle. 

We provide in this section the proof of the corollaries that showcase appropriate parameter- 
tunings for achieving sublinear query time and subquadratic preprocessing requirements. 


Proof of Corollary Again we set w = in order to achieve (sublinear) query time 
n^. We then set ^ = 1 + a ■ {1 — 0) ^ 0 = 2)^+0 ■ Observe that, for this value of 9, it 
is guaranteed that | < 3, as was assumed in the proof of Lemma To guarantee 

subquadratic preprocessing requirements n '^~we must assure that: 


<tP- 


2 — /3>2 — w-PaAl — 0)<t^/3<w — a-(l — 0)<^/3<- a ■ — - 

r + 1 2/v + a 

r < I ■ T_ 1 

- “ l(S+“)+(S-i) 


The approximation guarantee is, again, the one provided by the RQA query algorithm. ◄ 


Proof of Corollary |11[ Recall that, in order to assure a stretch factor 1 + k ■ e, we must set 


the recursion budget r > 


log log(l+£/'l/)) 


— 1 = ri{k), as in the proof of Corollary 


From Theorem we also have an upper bound on the recursion budget. Thus, as /3 0, 

a sufficient condition for the recursion budget, so that the required stretch is achieved is the 
following: 


r 7 (/c) 


2 <^ 

a 


2lv- 


a 


2lv- 1 


a < 


25 


— [ri{k)+2]-{2—i^) — 5L' 


It is now straightforward that the expected query-time is indeed whereas the pre¬ 
processing requirements are since we consider a very small value for (3. ◄ 


F I Analysis of the HORN oracle 

The construction of the travel-time summaries for HORN is based on the FLAT (BIS -P TRAP) 
preprocessing scenario. The queries are served by the HQA query algorithm. The oracle 
exploits two fundamental properties: 

(i) the approximation guarantee of a path via some landmark £ strongly depends on the 
relative distance of the landmark from the origin o, compared to the distance of the 
destination d from o; 

(ii) given that the expected distance of a level-z from the origin is roughly jt is rather 
unlikely that the first level-* landmark will appear too early or too late (i.e., outside a 
sufficiently wide ring-stripe around (o, to))- 
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Property (i) is exploited by the (ESC) criterion in order to handle the exceptional case 
where a higher-level landmark (which also happens to be informed) appears before the first 
informed landmark from the appropriate level. Property (ii) is actually an event that holds 
with high probability, as is shown in the detailed analysis of the oracle, and is exploited 
by the (ALH) criterion. Therefore, the event that an informed landmark appears which 
is also at the right distance, whereas the previously discovered landmarks (most likely of 
smaller levels) were uninformed, reveals an asymptotic bound for the unknown Dijkstra- 
Rank r[o, d]{to) of the destination. 

We now provide a sequence of lemmata which will eventually be used in the proof our 
main technical result, concerning the complexities of HORN mentioned in Theorem 12 (cf. 
Section . 

We start with an upper bound on the free-flow distance of a discovered landmark £o from 
d. At this point we do not require that d € C[£o\- 

► Lemma 15. Let G L be a landmark diseovered by HQA. Then it holds that i2[^ojd] < 
ini— ■ D[o,£o]{to) + D[o,d\{to). 

Proof of Lemma 1151 By Assumption |2.2| we know that: 

D[lo,0]{to)<C-D[0,lo\{to) (1) 

By Assumption |2. 1 1 we also know that: 


Va; > 0, -Ami„-x< D[to,o\{to) - D[io,o]{to - x) 

^ D[£o,o]{to-x) < D[to,o\{to) + h. min ‘ ^ 


( 2 ) 


We look for a particular departure-time to — Xo, and the corresponding minimum-travel-time 
D[£o, o]{to — Xo), so as to be at the origin o exactly at time to- That is: 

to - to Xo “t" L)^£o^ ^](^0 Xq^ 

/. 0 ./ 

Xo — L)^£o^ o](to ^o) ^ L)^£o-) o](to) -f A,^in * Xo 




l-Ao 


1- An 


(3) 


Finally, we upper-bound the free-flow distance of £o from d by exploiting the triangle in¬ 
equality: 


D[£o,d] 


< D[£o,d]{to - Xo) < D[£o,o\{to - Xo) + D[o,d\{to) 

/• a */ /■ 

< -^- ■ D[o,£o]ito) + D[o,d]{to) 


1 - A„„„ 

which is exactly the desired inequality. 


The following lemma provides an upper bound on the approximation guarantee of HQA, 
when an (ifS'C)-termination occurs. 

► Lemma 16. Assume that io & L is informed (i.e., d S C[£o\) and settled by the initial TDD 
ball that HQA grows from {o,to)- Then, for any given value ^^>1, if HQA terminated due to 
oceurrenee of (ESC), the reported travel-time is at most a -|- e -|- -approximation 

of the minimum travel-time Rd = D[o,d\{to)- 
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Proof of Lemma [TEJ Let D[o,^o]{to) = Ro- Then, Rd < A[o,d]{to) = Ro+A[io,d]{to+Ro)- 


By Theorem 14 (cf. also the analysis of FCA in [TB]) we can easily deduce that: 


/* (£SC)-termination */ 


A[io,d]{to + Ro) ^ 

Ro Ro 

/I I \ / I i\ I ; 1 ^ A[£o, d](to +-Ro) ^ , ~^Rd , , 

--—--— < {1 + £)-— + ip -1 

Ro Ro 


Rd 

Ro 


>(p-{r+\) 


4’Ro 

l + £+ ^ < l + e + 
Rd 




(p ■ {r +1) 


Since 1 -I- £ -I- is an upper bound on the approximation guarantee provided by FCA (cf. 


Theorem 141, which is indeed simulated by HQ A until the determination of the appropriate 


level in the hierarchy, we conclude that the eventual solution that will be provided by HQA 
is at least as good, since the (ifS'C')-termination returns the best approximate solution seen 
so far via an informed landmark, among which is also the one that goes via to- 


We proceed now by studying the first appearance of a level-z landmark within the unique 
outgoing ball from (o,to)- The next lemma shows that, with high probability, this first 
appearance of a level-i landmark will take place in the following ring for level-i: 


RING[o-,i]{t,) := B [o ; • ln(n)l (to) \ R 


(7V.)V(>-+i) 

ln(n) 


{to) 


= B 


o ; 


ln(n) 
P* J 


(to) \ B 


o ; 


Pi ln(n)_ 


{to) 


since Ni = and pi = n i)/[(’'+i)7 ] = <5/(r+i)^ 


► Lemma 17. Vt G [k], there is at least one level-i landmark in RING[o]i\(to), with prob¬ 
ability 1 — o{^). 

Proof of Lemma ITTl Consider any subset of vertices S' C 1/, of size s = |S| € N. The 
probability that none of the vertices in S is a level-t landmark (i.e., from Li) is (1 — pi)® < 
exp(-spi). 

Observe now that, for i G [k], s, = |i?/A^G[o; t](to)| = Thus, we conclude 

that: ¥{\RING[o-,i]{to)nL,\ = 0}<exp{-s,-p,) = "^^^^^^GO{^) . ◄ 


The next lemma states that, given the actual Dijkstra-Rank r[o, (t](to) that we seek for, the 
level-t landmark (if any) that is settled within RING[o;i]{to) is indeed informed about d, 
for all i such that Ni > r[o,(i](to)- 

—1 

► Lemma 18. For i G [k], let r[o, d](to) < Ni = n . Assuming that {ALFI)-termination 
occurred, the first level-i landmark ii^o G LinRING[o;i]{to) that is settled by the initialTDD 
ball grown by HQA, has d G G[li^o\- 

Proof of Lemma [THl Let Ri^o = D[o,£i^o]{to) and recall that Rd = D[o, d\{to)- Assume also 
that Rd > Ri^o, because otherwise an exact solution will be anyway discovered before ii^o is 
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settled. Then we have, by Lemma [TS 

c 


D[ko, d] 


m,c,d] 


< 


< 

/. As. flU V 
< 

/. As. flU v 
< 


1-A„ 
1 + 

1 + 


Ri,o “t" Rd 


c 

c 

1 Amir 


' Rd 

■5(n) 


fin) ■ iD[e^,o, d])^ < f{n) • g^(n) • 1 + 


C 


1 ^min 

(1 + A)0( ^“fogT^rO . 77,'^ 1-A^in ) ! 

< n 


• n 


(7'-l)/7' 


n' ' \ / • n 

^(7'-l)/7'+c(l) < 


for any > (^+^y^°siog(u)+j.^iog{i+ certainly true since A G and 

log lQg(n.) \ 

f{n),g{n) G polylog(n) = n V A The last inequality implies that d G (^[A.o]- 


The next lemma states that at no level-j in the hierarchy earlier than the appropriate level 
i* corresponding to r[o, d](to) (cf. Theorem 12 1 , may HQA find a landmark which contains 
d in its coverage (and which would then cause an incorrect “guess” of the appropriate level 
for the query at hand), provided that no (ifS'C)-termination occurred. 


► Lemma 19. For i> 1, assume that r[o, d](to) > Ni = and, while growing the 

ball from {o,to), no (ESC) occurs. Then, VI < j < i no level-j landmark i in RING[o] j]{to) 
contains d in its coverage: d ^ Ui<j<i UfeLjnfi/AfG[o;i](to) C'[£]. 

Proof of Lemma 1191 Recall again that Rj o = D\o, (.j^o]ito) is the distance from the origin o 
to the first level-j landmark that we meet in RING[o] j]{to)■ We start by providing a simple 
proof for the cases of 1 < j < i — 1. We shall then handle the case j = i separately, since it 
is a little bit more involved. So, fix arbitrary j G [i — 1]. For any ^j^o G RING[o-, j]{to) H Lj 
it holds, by the triangle inequality, that: 


d](to + Rj,o) 


> 

/. As. flU v 
> 

> 


D[o,d\{to) - D[o,ij^o]ito) 

(r[o,d](to))^/^ -g(n)(r[o,4o](to)) 


1/A 


1 


(7^-1)/(A7A . 


,(7’-i)/(A7b _ gi^n) lnA^(n)n“(T''-i)AV+i)A7^j 


^( 7 ’ ^- 1 )/(A 7 ’) 


-gin)- 


ln^/^(n) 

,(l-FfT)(A-l)/(A7^') 


^ T[£j^o, d]{toRj^o) ^ 


/.As^^./ /D[£,,o,d]ito + Rj,o) 


\ gin) 

„(7-^-i)/(A7') ln^/^{n) 


-\ A 


_ g^in)f^C{n) ^a-1('^)„(i-4t)(7^-i)/(A7L_ 
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Observe now that, from this last inequality and Assumption |2.3[ the following is deduced: 

-\ A 


Uko.d] > 


n 


(7^-l)/(70 


polylog(n) 


n 


(7* J-1)/(A7’) 






> n 


-1^-1 , J-1 rnf \2 loglog(7i)\ 

-;3- + —i -■ log(r.) ) 


The last inequality holds when A G ^ ~ 7~* = The last 

inequality implies also that d ^ C[£j^o\- 

We shall now study separately the case j = i. Apart from the containment of the level-i 
landmarks in RING[o; i](to), we also exploit the fact that the first of these landmarks met by 
the unique TDD ball from (o, to) does not cause an (AS'C')-termination of HQA. This implies 
that: 


£^[£i,o, d](to + 


Ri. 


Ri 


Ri 


< 


> 


(1 + e)l3{r + 1 ) + '0 - 1 =: X 

£^[£i,oid\{to Ri^o) ^ D[ii^o,d\(to Ri^o) 


/* triangle ineq. */ 
> 


> 

A Ab. [T 4l «/ 
> 


X 

Rd Ri,o 
X 
Rd 

1 + x 

^(7‘-l)/(>7’) 


(i + x)/^/^H 

Nevertheless, we also know that for any G RING\o] i]{to) the following holds: 
r[o,£.,o](to) < -InH 

Ri,o = D[oJi^o]{to) < (r[o,£i,o](to))^/^ 

< c,(n)n‘5(T-i)/(7+i)A7b iniA(„) 


(4) 


A As. [^ ./ 


< 


7j(7’-1)/(>‘7’) 

(i + x)/^/^W 


(5) 


Inequality © holds if and only if 

7 * - 1 


1 - 


r + 1 


„(i-4T)(7‘-i)/A7b > {l + x)g{.n)f/^{n)W/^{n) 

log(n)-log(ln(n))-log(/(n)) > Alog(l + x) + Alog( 5 f(n)) 


1 - 


r + 1 


7 * - 1 


log(n) - ©(log log(n)) > A [log(l + x) + 0(log log(n))] 


which is certainly true for A G and f{n),g{n) G polylog(n). From (41 and 

© we are led to a contradiction. Therefore, any level-i landmark in RING[o]i\{to) either 
does not possess a travel-time summary for d, or causes an early-stopping of HQA upon its 
settlement. ◄ 


The last lemma proves that, given that an (ALiA)-termination occurs, HQA achieves the 
desired approximation guarantee. 
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► Lemma 20. Fori* € {2,3,..., fc + 1}, let < r[o,(f|(io) < Ni». Assume also 

that an {ALH)-termination oecurred. Then, an ^1 + —approximate solu¬ 
tion is returned hy HQA, with probability 1 — The expectation of the query-time is 

C>((Afj.)“polylog(n)). 


Proof of Lemma 1201 As it was already explained by Lemmataandwe know that 
a successful “guess” of i* will occur with probability 1 — 0(^)- Then, HQA proceeds with 
the execution of RQA^., whose expected time is 0((A^i.)“polylog(n)) (cf. Theorem 14l. 

As for the approximation guarantee, since the analysis for RQAj. is based solely on the 
quality of the paths via landmarks discovered from ball centers located along a minimum- 
travel-time od-path p* S SP[o,d\{to) (cf. corresponding proof for RQA in [I^), it suffices to 
prove that all the level-(> i*) landmarks discovered from ball centers which reside at the 
(unknown) shortest od-path, are informed about d. 

For an arbitrary ball center Wi*j G p* and its closest level-i landmark for some i > i*, 
^i,j G Li, let Tti* j — T)\,Wi* -\- D [o, (to)) ■ Then, either Ttri if Lii*,oF ■ * ■ T 

in which case an exact solution is returned, or else the following holds: 


D[^i,j,d] < 

A Ab. 
< 

A Ab. 

=>T\iij,d] < 

< 




+ Rd < ( 1 -|- 




Rd 


\ -L -^^min / 

f{n) ■ {D[£ij,d])^ < f{n) ■ fl+ ) • g^(n) • 

\ -L -i^min / 




for Alog( 5 (n)) -k log(/(n)) < C* log(n) <tk> > (A + 1) • Again, a sufficient 

condition for this is that A G o^^ ■ We therefore conclude that d G C[iij]. 

Since all the discovered landmarks from ball centers along the unknown shortest path 
are informed, the analysis of RQA,;. provides the claimed approximation guarantee. ◄ 


Having proved all the required building blocks, we are now ready to provide the proof of 
our main argument for HORN. 

Proof of Theorem 1121 We start with a short sketch of the proof, before delving into the 
details: The claimed preprocessing space and time are proved, mainly due to the appropriate 
choice of the fi parameters. We thus focus on the more involved analysis of the time- 
complexity of HQA. First of all, exploiting Assumptions |2.1| and |2.2| we prove an upper bound 
on the free-flow travel-time of the destination d from an arbitrary landmark £, as a function 
of the (unknown) travel-time iA[o,d](to)- We then prove that if the (AS'C')-termination 
criterion occurs, then we already have a satisfactory approximation guarantee. Otherwise, 
assuming that the Dijkstra Rank F[o, d](to) lies within the two consecutive values A^a-i and 
IVa, we exploit the triangle inequality along with Assumptions |2.3| and |2.4| to prove that, 
with high probability, the closest level-z* landmark £a,o to the origin o lies in a certain ring 
(i.e., difference of two cocentric balls from o) and is indeed “informed”. Additionally, all 
the landmarks discovered prematurely with respect to the ring of the corresponding level 
to which they belong, cannot be “informed”, given the failure of the (AS'C')-termination 
criterion. Moreover, discovered landmarks of higher levers, despite being informed, do not lie 
in the ring of the corresponding levels and therefore do not interrupt the guessing procedure. 
This then implies that the appropriate variant of RQA will be successfully executed and will 
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return a travel-time estimation with the required approximation guarantee and query-time 
complexity sublinear in Ni* (the targeted Dijkstra-Rank). 

We now proceed with the detailed presentation of the proof. For each level i G [fc], the 
targeted Dijkstra-Rank is W = . The coverage of each level-i landmark contains 

Ci = Ni ■ rfi' = +^i destinations, Fi = of which are discovered by BIS and 

the remaining Ci — Fi destinations are discovered by TRAP. 

We start with the analysis of the approximation guarantee and the query-time com¬ 
plexity for HQA. Observe that the query-time is dominated by the scenario in which there 
is an (AL7A)-termination (the discovery of the destination, or an (FFCj-termination will 
only improve the performance of the algorithm). Assume that for some i* € {2, 3,..., k + 


1} (which we call the appropriate _Jevel for the query at hand), = n 

2 _ 

T[o^d]{to) < n T*"' = Ni*. Lemma 


--^-1 

*-i 


< 


l + e + 




16 


assures that, when an (F5'C')-termination occurs, an 
^ (^(^ 1 ) j-approximate solution was anyway discovered and thus there is no need to 

£(r-i-i)^i_^ , the approx- 
Lemma 


guarantee the success of the guessing procedure. By setting ip = 
imation guarantee when (FS'C)-termination occurs is equal to 1 + £ ^/L/^y '+^-1 
proves that, with probability l — at least one appearance of a level-z* landmark occurs 

in the appropriate ring: Li* n RING[o; i*]{to) yf 0. Lemma 18 then proves that, given that 
an (ALiL)-termination occurs, the first level-i* landmark £i*^o settled by the initial TDD ball 
within RING[o; i*]{to) is indeed an “informed” landmark: d € G[ii* 


Lemma 19 


proves 

that, given an (ALiA)-termination, no landmark of a previous level j < i* that was settled be¬ 
fore £i*^o may be informed: V£ G Uj^[i*-i]{Ljr]RING[o;j]{to)), d ^ G[P\. As for landmarks 
I G if we settle i before and it happens that d G G[£] (i.e., £ is “informed”), 

this event will not interrupt the guessing procedure of HQA, because it is not in the ring of the 
corresponding level: (i?[o](to) j](to)) H j](to)) = 0. There¬ 
fore we conclude that, in the case of an (ALiL)-termination, with probability 1 — the 

“guess” of the appropriate level-i* is indeed correct. Finally, Lemma demonstrates that. 


{l+e/i^r+^ 

(l+e/'0)''+'--l 


when an (ALiL)-termination occurs with a successful “guess” of z*, a (1 -T e 
approximate solution is returned by HQA. 

As for the expected query-time of HQA, we wish to assure that it is sublinear in the 
Dijkstra-Rank of the appropriate level z* of the hierarchy. Suppose that the exponent of 
sublinearity is (5 G (a, 1), i.e., the query-time is comparable to 0(^{Ni*Y'), where a G (0,1) 
is the exponent relating the period of the metric with the network size (T = n“). We 
focus on the highly-probable event that HQA, if it terminates due to the (ALH) criterion, 
makes a correct guess of z*. Conditioned on this event, the expected cost of RQAj. gives the 
appropriate value for the exponent uii of the landmark-sampling probability: 


r+1 


log ( ^ ) polylog(n) = n‘^dr+i)+di) = 


'' r+l 7’’ 


The contribution to the expected query-time of HQA of the unlikely event that the algorithm 
makes a wrong guess about the appropriate level of the query at hand is negligible, due to 
the quite small probability of this happening. 

We proceed next with the study of the required time and space for the BIS -f TRAP- 
based preprocessing of the HORN oracle. We wish to bound the preprocessing requirements 
with (fc -I- 1) • for a given /3 > 0. We shall make the appropriate choices of our 

tuning parameters so that in each level it holds that the preprocessing requirements are 
F, , P, < rz2-/3. 
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We begin with the determination of the requirements of level-i, for each i G [k]. For this 
level we generate \Li\ = landmarks, each of which possesses travel-time summaries 

for all the Ci destinations contained in its own coverage. Based on the analysis of the 
preprocessing analysis of FLAT (cf. Theorem]^, but now restricting ourselves within the 
coverage of each landmark, we know that the overall preprocessing requirements of level-i 
are bounded as follows: 


S, € CI(|L,|.(F.^polylog(F,) + T‘-»'c,)) 


c O 




+°( 1 ) a(l-6»i)-|-2f^+J. 

+ n T 




) 


provided that < 'j b We choose again 9i = 27 TTT’ ’'^hich in turn assures that 


Si , PiG „2-p;.-Ha(l-e.)-ho(l) 


We correlate our demand for subquadratic preprocessing with the recursion budget (and 
thus, the approximation guarantee) that we can achieve: 




2 — uji — di) ^2 — <2 

P <uj^-a-{l-e,)= ^ - a- 


A (2/ig+a)(l-7 ') 
a P-{2/{av) + l)+2/i^-l 


Therefore, the overall preprocessing requirements of HORN are ^horn , 
= 77 ,^-^+TSio^F = 772 -/ 3 + 0 ( 1 ) ^ g polylog(n). 


Pmm G (^ + 1 ) ■ 

◄ 












